Sensor Deployment for Network-like 
Environments 

Luca Greco* Matteo Gaeta, Benedetto Piccoli^ 
June 18, 2010 

Abstract 

This paper considers the problem of optimally deploying omnidirec- 
tional sensors, with potentially limited sensing radius, in a network-like 
environment. This model provides a compact and effective description of 
complex environments as well as a proper representation of road or river 
networks. We present a two-step procedure based on a discrete-time gra- 
dient ascent algorithm to find a local optimum for this problem. The first 
step performs a coarse optimization where sensors are allowed to move in 
the plane, to vary their sensing radius and to make use of a reduced model 
of the environment called collapsed network. It is made up of a finite dis- 
crete set of points, barycenters, produced by collapsing network edges. 
Sensors can be also clustered to reduce the complexity of this phase. The 
sensors' positions found in the first step are then projected on the network 
and used in the second finer optimization, where sensors are constrained 
to move only on the network. The second step can be performed on-line, 
in a distributed fashion, by sensors moving in the real environment, and 
can make use of the full network as well as of the collapsed one. The adop- 
tion of a less constrained initial optimization has the merit of reducing 
the negative impact of the presence of a large number of local optima. 

The effectiveness of the presented procedure is illustrated by a simu- 
lated deployment problem in an airport environment. 

1 Introduction 

Imagine a scenario where a toxic gas is spreading in an area or a building and 
safe paths have to be found to evacuate people. Or think of an airport environ- 
ment where people moving through rooms and corridors has to be surveilled in 
order to detect and avoid terroristic actions. Or consider the need of measuring 
environmental quantities, such as temperature or humidity, on wide areas to the 
aim of improving theoretical models or making more accurate weather forecast. 

There is a great number of situations that would greatly enjoy the use of 
network of sensors. Indeed, many of the previous tasks are difficult, or impossi- 
ble, to be accomplished by a single sensor. The employment of a large number 
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of sensors increases the robustness to sensor failure and communication disrup- 
tion and make spatially-distributed observations possible. If sensors are able to 
move, the number of tasks they can perform is still greater. 

Static and dynamic sensors' networks need to be deployed in the environ- 
ment, and the way this problem is solved can significantly affect the quality of 
service they have to provide. 

1.1 Static Deployment and Locational Optimization 

Sensors' deployment problems are strictly related to facilities location-allocation 
problems, which are the subject of the locational optimization discipline (tHj)- 

In locational optimization objective functions are used to describe the inter- 
actions between users and facilities and among them. Users may find facilities 
desirable, hence they would like to exert an attractive force to facilities, or unde- 
sirable and they would repel them. The attractive model can describe allocation 
problems of useful services or facilities such as mailboxes, hospitals, fire stations, 
malls, etc. (see 0). The repulsive one, instead, can be used to model problems 
where polluting or dangerous facilities (i.e. nuclear reactors, garbage dumps, 
etc.) are to be located far enough from urban conglomerations. An excellent 
survey on undesirable facility locations problems is given by [T^] (see also [3]). 
These operational research problems can be converted in sensors' deployment 
problems by considering sensors as facilities and points or areas, where events 
can happen or some quantities has to be measured, as users. 

Two well known problems, involving one facility only, are the classical Weber 
and the obnoxious facility location problems (see [101 18] for a recent heuristic 
solution). Three problems involving p facilities are the p-center, p- median and 
p-dispersion problems. Some recent results on the p-center problem are in |29| 
and [5] . The latter paper addresses also the p-dispersion problem. 

A classical p-median problem, close to the one considered in this paper, can 
be simply described as the one of finding the optimal location of p facilities by 
minimizing the average distance of the demand points to the nearest facility 
(see for a recent survey on heuristics methods to solve it). Close to the 
p-median problem is the multisource Weber problem, for which many heuristics 
exist ([2). A more general formulation of these problems can be found in [7] and 
[TT] . In [7] a dynamical (gradient descent) version of the Lloyd's algorithm [23] 
has been presented to find a local optimum for a generalized p-median problem. 
A deterministic annealing optimization algorithm to solve the classical version 
is reported in [25] . The aforementioned solutions to the p- median problem, as 
well as many solutions to p-facilities problems, are based on the construction of 
a Voronoi Tesselation ([TT1[22])- 

1.2 Dynamic Deployment and Distributed Solutions 

The use of moving, instead of static, sensor networks provide a great flexibility 
in solving sensing tasks, mainly when the environment is partially or completely 
unknown or is not directly accessible for safety reasons. In these cases, sensors 
are usually initially deployed randomly and hence need to move in order to ac- 
quire knowledge of the environment and to optimally re-deploy for their task. 
Furthermore, environments are usually not static and the network may experi- 
ence sensor failure or loss. In these situations the properties of adaptivity and 
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reconfigurability owned by a network of moving sensors turn out very useful. 

A general tendency in robotic networks is to have sensors (agents) endowed 
with the same computational and sensing capabilities. This choice increases the 
overall robustness of the network, but usually calls for distributed coordination 
algorithms. Having equal sensors, indeed, naturally leads to define optimization 
and coordination algorithms based on local observations and local decisions 
([131 US [H]). Many of the algorithms proposed in the previous section involve 
the solution of a global optimization problem requiring a complete knowledge 
of the environment and of sensors' distribution. The solutions to p-center, p- 
dispersion and p- median problems proposed in [71 [6] , instead, are all spatially 
distributed, with the meaning that each sensor requires only the knowledge of 
positions of its neighbors (or even less if it has a limited sensing radius) . This 
fact allows a distributed implementation where each sensor computes its next 
movement without centralized coordination. 

Other solutions to the area-coverage problem look at sensors like particles 
subject to virtual forces or potential fields. The compositions of suitably defined 
attractive and repulsive forces is then used to make the network behaving in 
the desired fashion (spread sensors, avoid obstacles, keep connectivity, etc.). 
Representative for this kind of approach are the algorithms presented in [T71I35] . 
or in [21], where also secure connectivity issues are considered. In [TB], instead, 
it is raised the relevant problem of power consumption in wireless networks and 
three energy-efficient algorithms are presented for sensors' deployment. 

1.3 Network-like Environments and Paper Contributions 

In this paper, we focus on network-like environments as there are surveillance 
or monitoring problems where such a kind of model can provide a more suitable 
description. 

Network models represent a natural choice whenever environments have an 
intrinsic network structure. It is the case, for instance, when sensors have to 
be optimally deployed over a network of roads to monitor vehicular traffic, or 
in a river network to measure temperature or pollutants concentration. Even 
some location-allocation problems can involve networks. Consider, for instance, 
the case where useful facilities (i.e. schools, hospitals) have to be located in 
the interior of a network of roads, which is the source of a nuisance (i.e. noise, 
pollutants) , with the goal of minimizing its harmful impact on them (see [B] for 
the case with one facility) ; or the dual problem of locating obnoxious facilities 
(i.e. dumps, industrial plants, mobile phone repeater antennas) reducing the 
hazard on the network. 

Most notably, we think that a network-like model can provide an effective 
and compact description for complex environments, focusing only on major 
features and abstracting from those geometrical details which are less important 
for the deployment problem. The coverage of nonconvex environments with 
holes or obstacles, for instance, is a challenging task (U), which can enjoy 
significantly the use of a reduced network-like model. Environments with a 
complex structure, accounting for a large number of variously sized and shaped 
rooms, passages, forbidden areas and obstacles, can be reduced to a set of 
connected paths where the sensing task is more requested or where sensors are 
forced to pass. An airport is a very representative example of this kind of 
environments. In this case people moving throughout the airport can be aptly 
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compared with a network flow and focus can be on paths more than on corridors, 
hahs and lounges. 

Many of the problems introduced in previous sections have been formulated 
even for a network- like environment (^ [301 131j). but they usually consider a fi- 
nite discrete set of demand points located on the network's nodes and try to 
optimize the locations of facilities w.r.t. some objective function accounting for 
the distance from them. An important thread of works for the deployment of 
sensors both on plane and on a network is represented by the papers [121 [101 [21] ■ 
In such works, the authors propose powerful greedy algorithms that provide a 
constant approximation of the optimal solution. Their method, however, aims 
at solving a global static deployment problem, considering a finite set of demand 
points and allowing sensors to jump among positions. 

In this paper, instead, we consider a deployment problem where sensors 
use local information to dynamically solve the optimization problem while they 
are moving in the environment. The constraints induced by the environment 
to the motion of sensors are explicitly considered. More precisely, we address a 
generalized p- median problem involving omnidirectional sensors with potentially 
limited sensing radius and we extend the formulation presented in |S1 [71 [5S] to 
network- like environments. The task is to find sensors' positions that optimize 
an objective function defined on the network and accounting for the sensor's 
features and preferential areas. This is a mixed problem, since the network is 
considered embedded in the plane (it is a continuous set of demand points) and 
the planar euclidean norm is used to measure the distance between sensors and 
network. 

The core of the cited formulation and of our solution is a discrete-time gra- 
dient ascent algorithm based on Voronoi partitions and aiming at maximizing 
the objective function. It is a well known fact, however, that such a kind of al- 
gorithms can get stuck early in local optima, especially when sensors are forced 
to move in an over-constrained environment like a network. Moreover, the local 
maximum found by the algorithm is often greatly related to the initial sensors' 
position. 

For these reasons, we present a novel two-step procedure performing an 
initial coarse optimization, whose purpose is to provide a good starting point for 
a second finer optimization. The first step can be carried out off-line, either by a 
central unit, or by each sensor individually (without doing real movements). The 
impact of local optima is reduced by allowing sensors, in the initial optimization, 
to virtually move in and to vary their sensing radius arbitrarily. In order 
to reduce the complexity of this phase, sensors can be initially clustered and 
the optimization problem solved for the clusters' centers. After that, a desired 
number of sensors is spread close to clusters' centers and the sensors' positions 
thus found are projected on the closest edges of the network. The projected 
positions are then used in the second optimization, where sensors are constrained 
to move only on the network. The second step can be performed on-line, in a 
distributed fashion, by sensors moving in the real environment. 

The use of a two-step procedure is motivated also by the very nature of some 
surveillance tasks, such as for instance airport surveillance, where a large num- 
ber of individuals (sometimes referred to as mass objects |15| ) are monitored. 
In these cases, sensors can solve the first step optimization using imprecise or 
estimated information and keeping still; then they can move to reach the final 
projected positions using planned routes compatible with the network. After 
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this initial deployment, sensors can change their positions by dynamically solv- 
ing a distributed optimization problem (second step) based on real measures 
taken from the (potentially varying) environment. 

It is worth noting that the present procedure can be used to solve both 
static and dynamic deployment problems. Moreover, the first step deserves 
attention by its own, since it provides a solution to those problems involving 
facilities located in the interior of the network mentioned at the beginning of 
this section. 

Another contribution of this paper is the introduction of a simplified model 
of the network (similar to the discretization in [14] ) called collapsed network and 
consisting of finite many points. It is obtained by decomposing each segment 
of the original network in one or more sub-segments and collapsing each sub- 
segment in its barycenter. This model allows a coarser but faster optimization, 
since computations with barycenters are remarkably less than those needed by 
the full network. Collapsed network, hence, is intended mainly for fastening the 
first optimization, but can be used profitably also for the second step. Indeed, 
it turns out particularly useful in practical implementation involving hardware 
with limited computational capabilities. 

As mentioned above, our work is related to that of [HI [Zl [2S]. In particular, 
the first step of the optimization, allowing sensors to move in M^, could be re- 
garded as a specialization of the problem described in [7]. However, the different 
topology induced by the network introduces issues related to the explicit com- 
putation of the gradient and to the convergence of the maximization algorithm, 
which deserve special solutions. A relevant difference is that the gradient of the 
objective function presents discontinuity points caused by barycenters on the 
boundary edge of two neighboring Voronoi cells. Such barycenters can change 
allocation during sensors' motion, inducing abrupt variations to the value of the 
objective function associated to each cell. This fact prevents a classical con- 
vergence proof for gradient algorithms, hence we consider our proof as a minor 
contribution of the paper. Some results about convergence may alternatively 
be derived by using the method of Kushner and Borkar ([33J [T]) of stochastic 
approximation to deal with our differential inclusion. 

The outline of the paper is as follows. In section [5] the mathematical defi- 
nitions of sensors, network and Voronoi covering are introduced along with the 
objective function to be maximized to solve the deployment problem. Section [3] 
is devoted to the introduction of the collapsed network and to the formulation, 
and proof of convergence, of a gradient ascent algorithm to solve the first step 
fsubsection l3.ip and the second step (subsection [321) of the optimization proce- 
dure. Section m addresses the solution of the second step optimization involving 
the full network. Finally, a network describing an airport environment is used 
in section [5] to illustrate by simulations the effectiveness of the proposed opti- 
mization procedure. Conclusions and future research directions are reported in 
section [6] 

2 Preliminaries and Problem Formulation 

In this section we introduce the mathematical framework to describe the sensors, 
the network and its Voronoi covering. 

Definition 1 Given two points pi,P2 G I^^j with pi ^ p2, si2 = [pi,P2] C 
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is the segment joining pi and p2 and — (pitP2) is the open segment between 
them. We define length of a segment S12 as l{si2) — \\p2 —Pi\\, where \\-\\ is the 
Euclidean norm; barycenter of a segment S12 the point b{si2) = \ {p2 +Pi) S 
S12; partition of a segment s = [pi,p2] in k sub-segments, the set of segments 
{siL=i,....fc 9^ven by 

\ . f. -,^ (P2-P1) , . {P2-P1) ' 
Si^ Pi + {i~l) ,Pi+i ^ ■ 

Definition 2 A network M = (V, S) is a subset of consisting of a set of 
points V = {wi, . . . , u„ 6 R^, ^ Vi 7^ j} and a set of segments S C {s,j — 
[vi,Vj\ C M^, i, J G {1, . . . ,ri} i ^ j}, such that: 

i) Vui e V, 3vj G V , Vi ^ Vj such that Sij d S (no isolated vertex); 

ii) V«, j, h,k ^ {1, . . . , ri} , (i, j) ^ (/i, fc), s^jHs^j;, — 9 (no segment intersection) . 

Definition 3 Given a network Af and a set of points V — {pi, . . . ,Pm} C TV, 
the Voronoi covering of J\f generated by V with respect to the Euclidean norm 
is the collection of sets m} '^'^/^^^c? by 

Vt^iV) = {qeN\\\q-p.\\< \\q-p,\\ , Vp, £ V} . 

Remark 4 It is straightforward to recognize that V^^iV) can be equivalently 
defined as V/^{V) = Vi{V)riJ\f , where Vi{T') is the i-th cell of the usual Voronoi 
partition o/R^ generated by V . The previous definition is about a covering and 
not a partition since neighboring cells can have a nontrivial intersection: a 
portion of a segment can belong to the shared edge of two cells Vi('P) andVj{'P). 

We adapt the framework provided in [B] to describe the sensors' and network 
features. Each sensor is modeled by the (same) performance function f : K+ — > 
R, that is a non-increasing and piecewise differentiable map having a finite 
number of bounded discontinuities a.t Ri, . . . , Rn £ R+, with i?i < . . . < Rpf. 
We can set Rq = 0, Rn+i — +00 and write 

N+l 

fix) = J2 /a(a:)l[i^_„fl,„)(x), (1) 

Q = l 

with fa : [Ra-i, Ra] ^ R, « G {!,..., 1} non-increasing continuously 
differentiable functions such that fa{Ra) > fa+i{Ra) for a G {1, . . . , N}. In 
order to model regions of the network with different importance, we can use a 
density function (j) : J\f ^ R+ , which is bounded and measurable on TV. Given 

(7 : R^ — i> R we indicate by / g{q)dq (respectively / g{q)dq) the sum of the 

J^f Jv^^ {V) 

linear integrals of g over the segments of TV (respectively {V)) using an arc- 
length parameterization. With these functions we can define the multi-center 
function % : TV™ — !> R for m sensors located in V = {pi, . . . ,Pm} C TV 

H{P)^ f raa^ fi\\q^p,\\)4>iq)dq. (2) 
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We can also provide an alternative expression for ^ based on the Voronoi 
covering induced by V as follows 



n{r)^Y. nh-p^\\)'P{q)dq- / nh-Ph\\)'P{q)dqAi) 



where Af^ ^ V,^ {V) n V^{V) and ^ {Af^ | h < fc, V/i, fc G {1, . . . , m}}. 
The second term in ([3]) is not null if and only if there exists a non trivial segment 
s C Sij G S such that s C A-^, for some i, j G {1, . . . , n} and ft., G {1, . . . , m}. 



3 Deployment over a Collapsed Network 

In a collapsed network each segment of the original network is decomposed in 
one or more sub-segments and each sub-segment is collapsed in its barycenter. 
Chosen a value for r guaranteeing a good approximation, we can build the 
r-coUapsed network as follows: 

Definition 5 (r-Collapsed Network) Given a network N = (V,5) and r > 



0, Vs G consider its partition in kg 



e{s) 



sub-segments Si (having at most 



length r) and the associated set of barycenters {b {si)}^^-^ ^ . We define the r- 
collapsed network associated to Af the set of points = U (■^Oli^i 



The multi-center function must be re-defined since the integration domain 
is now a discrete set represented by the barycenters. Hence we set 

niV)= max (4) 

where (pb^ are suitable (density) weights assigned to barycenters. 
3.1 Sensors Moving in 

In this section we solve the deployment problem for a collapsed network and 
sensors moving in M^, which constitutes the first step of our optimization pro- 
cedure. 

Also for the multi-center function Q we can provide an alternative expres- 
sion using the Voronoi covering. We need the following dcfinitiorQ: 

Definition 6 Given an r- collapsed network for some r G and a set 
of points V — {pi,...,pm} C M'^, the Voronoi covering of generated by V 

with respect to the Euclidean norm is the collection of sets \ V^'' [V) \ 

L J ie{l,...,m} 

defined by 

vf^iP) = {beCf\ \\b-p,\\ < \\b~p,\\ , Wp, G V} . 



^For this definition similar remarks as Remark 2] also apply. 
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We define also the boundary of a Voronoi cell as 

dV^^iV) ^{be Cf I II& -p.ll = Wb^pA , Vp, e V} , 
and, in order to simplify the problem, we make the following assumption 

Assumption 7 dV^^ {V) = Vi £ {1, . . . , m}. 

With this assumption the multi-center function ^ can be written also as 



(5) 



b^evf'' {V) 



or, setting dist((7,P) = minig{i^...^„i} Hg as 



(6) 



Remark 8 'H{V) is not globally Lipschitz as /(•) is not a continuous func- 
tion. However, if f{x) is continuous and piecewise differentiable with bounded 
derivative, then H('P) is globally Lipschitz. In order to prove the global Lips- 
chitz continuity, let us consider two sets of points V — {pi, . . . ,Pm} C M'^ and 
V ~ {p'n . . . tP'„,} C and compute 

n{V)-H{P')^ J2 [/(dist(6e,P))-/(dist(6e,7''))]0f,e 

df 



^ E 



dx 



< 



where 



^b^ec^ 
is the Loo-norm of jJ-. 



|dist(6e,7') -dist(6e,7'')l'/'fce 
df 



dx 



\v-r'\\, 



Theorem 9 The multi- center function H is continuously differentiable on 
{T)qmY^ , where 

^c^^ - U^ec^ i'? ^ I 11^- - ^ll = ^" Vz = 1, . . . , TV} 
is the discontinuity set of /(•) in R^. Moreover, for each h S {1, . . . , m} 



dniv) 



d 



dph ^ dph 



E n-/(ll&e-p.||)0.. 



(7) 



Proof. The continuous differentiability of % on (R^)™ \ (Dcat)'" is a straight 
consequence of the same property of /(•) on 'Dqm . As concerns the gradient. 



8 



using ^ we have 



= E /-/(dist(6e,7'))06,. 

In the hypothesis of Assuniption[71 we have that the index set /f,e — argmin^gj]^ 
has cardinality 1. Hence, V6e such that I^e ^ i =^ h, -^f (dist(6e, P)) = and 

dH{V) ^ d 

^ {b,ecM/6e=/i} 

whereby the thesis follows using the definition of V^^'' {V). ■ 

The sensors' location-allocation problem can be addressed by means of a 
gradient-like algorithm. If a continuous time implementation is looked for, the 
following fictitious dynamics would be associated to the sensors' positions 

V = vniv). (8) 

Unfortunately, this dynamics conveys some problems. It is well defined as long 
as Assumption [7] and Theorem [5] are fulfilled, but these hypotheses are, in fact, 
too stringent for the algorithm to work properly. Indeed, they would require the 
evolution of the sensors to avoid any position in the discontinuity set and the 
barycenters not to enter or exit the Voronoi cells where they are at the initial 
time instant. 

First of all we can reduce the analysis to continuously differentiable functions 
to avoid issues related to the existence of the gradient. Still the relaxation of As- 
sumption[7]induces some problems on the definition of the gradient. Barycenters 
on a boundary edge of a Voronoi cell belong to all the cells sharing that edge. 
All sensors' positions producing these configurations are discontinuity points 
for S/H{V). Roughly speaking, the gradient takes different values depending on 
which cell the shared barycenters are assumed to belong to. This fact makes the 
equation ([8]) a set of differential equations with discontinuous right-hand side. 

The problem of shared barycenters can be solved by adding a lexicographic 
criterion to the definitions based on the euclidean distance. Indeed, with this 
criterion barycenters on boundary edges are allocated univocally to the sensor 
having the lower index (w.r.t. the Lexicographic Order (L.O.)) among sharing 
sensors. This fact allows us to define a genuine Voronoi partition, no longer a 
covering, whose generic cell is given by (compare with Definition [5]): 

= {6 e Cf I \\b-p,,\\ < \\b-p,\\ Vp, eP A 

\\b - p,\\ < \\b - pjW if j < i w.r.t. the L.O.}. (9) 

We must now define a generalized (lexicographic) gradient of H, WiH{V), 
according to this new definition. V6e G V^'- {V) \ dV^ {V) we use the classical 
formula given by ([7]). Vfee G dV^ {V) notice that the partial derivative of /, 
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'^;^f iW^e — Ph\\), exists and is well defined. In the light of this remark we can 
write the ft.-th component of the generalized (lexicographic) gradient of H as 

fWp. E ^Jm-.SK^ (10, 

which is formally equal to the formula (O provided by Theorem 1^1 

The differential equation using this new definition for the gradient, however, 
does not imply the existence and uniqueness of the solution, and this proof may 
turn out to be complex due to special sensors' and barycenters' configurations. 
Moreover, the formula (jlOp accounts only for the infinitesimal perturbations of 
sensors' position not inducing any barycenters to enter or exit Voronoi cells, 
hence changing their allocation. 

In order to simplify the convergence proof and to provide an algorithm which 
is more suitable for a realistic implementation, we consider here a discrete-time 
version of the gradient algorithm. In our case, the discrete-time implementation 
can overcome the problems with discontinuity gradients due to the properties 
of the function H and its discontinuity points, as it is shown by the following 
theorem. 

Theorem 10 Consider the following discrete-time evolution for the sensors' 
positions 

■p{k+i) ^ ^{k) ^ g^^^^^^(k)^^ 

where the h-th component ofViH is given by 110^) and T-i : K^™ — ^ R as m 

If f{-) has locally bounded second derivatives, then, for suitable 5k, V^'^^ lies in 

a bounded set and 

i) 'H{V^''^) is monotonically nondecreasing; 

ii) ■p'-'"'^ converges to the set of critical points ofH. 

Proof. It is easy to see that there exists a ball B ^ Af such that, if pi e dB, then 
^'^^^-^ points inside B. Thus the fact that T'^'^^ is bounded for 6k sufficiently 
small, easily follows. According to the definition in equation ([3]) of a Voronoi 
cell using the lexicographic rule, we can define 



be&Vp- {V) 



hence we can write T-l{V) = '^^'Hp^{V). We must prove that Hp. ('P^'^"'"^-') > 

1=1 

'Hp^{'P^'^^) for \/pi G V''^^ obeying the discrete-time evolution ([TT|) and that 
3p,'e pC^) such that 'Hp.(p('=+i)) > Up^V^'"^) iiV'^^^ does not belong to the set 



be - P\ 



of critical points of U. Define (pC^+i)) = ^ f 
that is the cost Hp-i-) computed by using the new sensors' positions pC^+i) but 
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the old allocation of barycenters to sensors, i.e. the Voronoi partition generated 

by T'C^) . Therefore, we write 



dpi 



•(V;?^(pW)).j+o((4)2) 



(12) 



where the (generalized) gradient is computed by means of the lexicographic 
assignment based on the Voronoi partition generated by P^*^). If the Voronoi 
partitions generated by •p('=+^) and V^'^^ coincide (same barycenters' allocation), 
then ?^p,(P('=+i)) = ^^^.(Pe^+i)) and we can assert directly that 'Hp^{V^''+'^^) > 

^p.(p('=)) or •Hp,(p('=^+i)) > •Hp,(p('=)) if ^^"'qp^"^^ + 0, that is if V^^^ does 
not belong to the set of critical points of 7^. It is worth noting that similar 
remarks about the strict inequality apply also in what follows and they will 
not be repeated. If the barycenters' allocation is different, 'Kp^iJ^'^^^^^) may be 
smaller than V.^.i^'^^^^^) and we cannot say anything about its relation with 
^p.('p('=)). This fact is due to the presence of barycenters changing allocation 
during the; sensors' evolution, hence the cells involved in the change cannot be 
considered independently. Let us consider for simplicity only one barycenter 

6e e Vr^- (P^*')) and suppose that at step fc + 1 &e e Vi"^ (P^'^+i)). With this 
assumption no barycenters can enter or exit the union of the two cells but 6e- 
We have 



Hp,{V^^+^^) 



he 


(fc+1) 
-Pi 


be 


(fe+1) 
-p) 



(13) 



= ?^^,(P('=+1)) + /(| 

= :W^.(p('=+i))-/( 

In other words, 'Hp.{P^''~^^^) is grown w.r.t. the ideal value 'Hp.{P^''~^^^) due to 

the allocation of 6e topf whereas, Hp^ (-p('=+i)) is decreased w.r.t. ('P('=+^) ) 

by the contribution that be would have given if it were allocated to Pj*^^^^ 
even at the step fc + 1. It is worth noting that if be changes allocation from 



P 



(fc+i) 



to p. 



(fc+i) 



then 



■P. 



(fc+i) 



< 



be 



P 



(fe+1) 



(the equality holds only 



if be e dVi'- (T'C^+i)) n dVj^ (7?(fc+i)) and, being i < j, the allocation is in- 
duced by the lexicographic rule). Therefore, due to the monotonicity of /, 



(fe+1) 



Pi 

UpSV^''+^'^ 



)>-<{ 



■P 



(fe+1) 



^ . Summing up, we have 



+ Hp^ (pc^+i)) > n;^ (pc^+i)) + h;^ (7'('=^ 



np,{p<^''^) + si 



np,{v'^'^) + 5i 



There exist some constants c^-' ,5^^ > such that 



dUp^iV^^'^) 



dpi 
dHp.iV^^'^) 



dpj 



+oi{sm 



+ 



n 



^^(^(/c+D) + ^^^_(p(fe+i)) > 7^p^(77W) + Hj,^.(p(fc)) + 5^^^ + o . 
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Extending the same reasoning to more complex configurations involving more 
than one common boundary edge and more than two neighboring cells, we can 
say that 3ck, Sk > such that 

m 

^(p(fe+i)) ^ > n{v^''^) + SkCk + o{sl). 

This proves assertion i). Assertion ii) can be proved by exploiting the results of 
Section 3.2 of [B]. More precisely, using the fact that P*^*^^ is bounded and / has 
locally bounded second derivatives, then there exists ^ > such that we can 
choose Sk > S. Then we conclude, using Proposition 3.4 of [6 , that assertion ii) 
is true. ■ 

Remark 11 (Distributed Implementation) The use of a gradient ascent 
algorithm based on a Voronoi partition, allows us to solve not only a static 
deployment problem, but also a dynamic one. As shown in JBjl, this kind of 
algorithms is spatially distributed, with the meaning that the i-th sensor needs 
only to know the position of its neighbors in order to determine the boundary of 

O'hi 

its cell and, hence, to compute . For the same reason the i-th sensor can 
choose the value of the step-size 5\ independently of the other sensors simply 
performing locally a classical line search algorithm. This property makes the 
algorithm suitable for a spatially distributed implementation. 

The independence in the choice of the step-sizes 5\. is obviously preserved 
in each period k, as long as a synchronous implementation is considered. In 
this case sensors have access to a global clock, or perform a synchronization 
algorithm. At the beginning of the the k-th period (instant tk), all sensors are 
idle, build their Voronoi cells and compute their gradients and step-sizes, then 
they move until, at most, the end of the period (instant tk+i). If, instead, an 
asynchronous implementation is considered, further hypotheses are necessary 
to ensure that independence is preserved. Unfortunately, the discontinuity of 
the gradient prevent us from using the results of 132]. But, if a sensor has 
the capability to detect when its neighbors start and stop moving and when a 
new sensor joins the neighborhood, the asynchronous algorithm presented in ^ 
(Table IV), can be applied, thus automatically recovering the independence. 

Remark 12 In the previous theorem, for sake of simplicity, we did not consider 
degenerate configurations where different sensors have the same position (pi = pj 
for i ^ jV- caw be proved that if the initial positions of sensors are not 

degenerate, sensors can always choose a suitable to avoid the occurrence of 
these configurations. 

3.2 Sensors Moving on the Network 

This section is devoted to the case of sensors constrained to move on the network 
and sensing a collapsed network. Therefore, these results are suitable for an 
implementation of the second step of our procedure on hardware with limited 
computational capabilities. 

We still assume /(•) to be a continuously differentiable function and we make 
use of the lexicographic criterion for the barycenter allocation. As concerns 
sensors' motion, however, we cannot use directly the gradient since the sensors 
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have to remain on the network. We must consider now the directional derivative 
of % along the edges of the network. 

Following the guidelines of the previous section, the following theorem can 
be proved. 

Theorem 13 Given a network J\f — (y,S) and the related r-coUapsed network 
C!^ , the multi- center function Ti is continuously differentiable almost everywhere 
on A/""*. In particular, on each open segment s°j such that Sij G S, given the 
unit vector Wij such that Sij ■ Wij = \\sij\\, the directional derivative in ph G s°j 
along Wij is 

D^,^n{r)[pu] = (^(^) • ^^.) (14) 
E ^fi\\be-PH\\)K- 

dph ^ dph 

b.ev^^ (V) 

The directional derivative is a multivalued function on the vertices of the 
network as more than one edge can share the same vertex, but we need a univocal 
definition. Hence, we fix a choice rule such that the directional derivative in a 
vertex is given by the maximum among all the derivatives defined for each 
possible direction that does not lead the sensor out of the network. If all the 
directional derivatives in a vertex point outward the network, then the derivative 
is set equal to zero. 

Definition 14 Given the set 

Svi = {s G iS I 3vj G V, 3(5 > s.t. s = [vi,Vj] V s = [vj,Vi], 

yS G [0,6] V, + 5D^^^-H{V)[v,] G N] , (15) 

we define the directional derivative of 'H['P) in any point ph (z J\f as follows 

Dn,,M'P)[Ph] g^^en by (I^ Vp,, G \ V 



max D^^MV)[ph\ = e V 

5., ^0 (16) 

yph = G V 



We can now define the discrete-time gradient-like algorithm. 

Theorem 15 Gonsider the following discrete-time evolution for the sensors' 
positions 

^(k+i) ^ j,(k) ^ SkOniv'-''^), (17) 

where the h-th component of DH is given by ild)) and T-L : TV™ -^M. as in 1^. 
If f{-) has locally bounded second derivatives, then, for suitable 5k, P^'^^ lies in 
a bounded set and 

i) %{V'^^^) is monotonically nondecreasing; 
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ii) •p(fc) converges to the set of critical points ofH. 

Proof. The proof is essentially the same of the Theorem [TOj except for the use 
of the derivative ([T6)) in place of the gradient ([T0| . In particular equality (fT2)) 
becomes 

for Pi moving along Whi ■ ■ 

4 Deployment over a Full Network 

In this section we consider a more accurate version of the second step of the 
optimization procedure, namely sensors constrained on the network and sensing 
the full network. To start with, let us define the boundary of a Voronoi cell 
as dV^{V) = {g G A/" I II9 — Pill = ||g — P^ll , ^Pj G V}, and the instantaneous 
discontinuity set of /(•) as 

T^N{V)^[j ^{q^M\ ||g-p,||=i?„ V* = l,...,iV}. 

Assumption 16 We make the following assumptions: 

i) orthogonality assumption: \fh,k G {!,..., n}, \/i £ {l,...,m} and for any 

segment s = [a, b] C s^k G S with a ^h, s <f_ 8X^/^(7^); 

ii) dV^^iV) n V^{r) = 0, V^ G {1, . . . , m}; 

iii) V n V^f{r) - 0; 

iv) V/i, fc G {1, . . . , n}, V^ G {1, . . . ,m}, Wq G ShkCW^^iV), if {q-Pi)-{vh-Vk) = 

0^ ||9-p,;|| i {i?i,...,i?jv}; 

With the orthogonality assumption the expression ([3]) simplifies to 

m „ 

nv)^Y. f{h-pA\)^{<i)dq. (18) 

Theorem 17 Given a network Af ~ (V,iS) if Assumvtion \TE\ holds, the multi- 
center function T-L is continuously differentiahle almost everywhere on A/"™ . In 
particular, on each open segment s°j such that Sij G S, given the unit vector Wij 
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such that Sij ■ Wij — ||sjj||, the directional derivative in ph £ s°j along Wij is 



D^,^,n{V)[pH] = ( — (P) • m,, ) w., (19) 
dU 

Oph 



rW ^^^^ 



k=l 



Tu = \H-a,\\ f ^''fl <P{nit))dt 

Jlos]\{tl„tl„...,t%„t%,}di^ IK-7feWll 

f^i ptlh-7fe(ia,,)|| 

where jk(t) = + (6^: — ak)t, t G [0,1] is a parameterization for the k-th 
segment [ak,bk\ S Vy^{V), Mh{V) is the number of segments in Vj^{P) and 
t'l j e [0, 1], j e {1,2} are the zeros of \\-/kit) - Ph\\ - = (if any). 



Proof. Consider the gradient of ^{{V) in the form 

OPh "Ph~^Jvf{v) 
d 



f{\\q~Ph\\)<t>{q)dq 

OPhJv^(V) 

+ ^E/ /(ll9-Kll)0(9)d<Z. (20) 



Let us consider the second term of (PO]) for each i ^ h 
d 



. / f{\\q-p^\\)<l^{q)dq 

OPhJv^(V) 



lim^l/ f{h-p^\\)(^{q)dq- i f {\\q ~ P^\\) {q) dq 



= ]im -^(f f{\\q-p^\\)cf>{q)dq- f f {\\q - P^\\) cb (q) dq] , 

(21) 

where AV^{V)+ = {{pi, . . . ,ph + e, . . . ,pm}) \ Vt^{V) and 

[V)- = Vf{P) \ Vfiipi, . . . ,ph +£,... ,Pm})- It is worth noting that 
V^{{pi, . . . ,ph + £,... ,Pm}) can be different from V^^ {V) if and only if pi £ 
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A/g^ (p/i, T-")!! Let us consider now the first term of (|20p 

-^f f{\\q-Ph\\)<lyiq)dq 
OPhJv^iV) 



limol/ f ih - iPh + e)\\) <P {q) dq - f f {\\q ~ Ph\\) <l> {q) dq 

lim ^ / (/ (Ik - {Ph + - / (Ik - PhW)) <P {q) dq 



limAl/ f{\\q-Pk\\)cl>{q)dq- I f {\\q - Ph\\) <f> {q) dq 



= lim ^ / (/ (Ik - (p„ + £)||) - / (Ik - P/.II)) (?) dq 
"^-^o llell-yy^^lP) 

+ lim^ / ifi\\q-{pi, + e)\\)-fi\\q-pn\\))<l^iq)dq 

- lim T^f (/ (Ik - (p,. + e)||) - / (Ik - PhW)) ct> [q) dq 

+ lim^(/ f{\\q-Ph\\)^{q)dq- [ f {\\q - Ph\\) 4> [q] dq] , 

(22) 

wfiere AVf^^{V)+ = F,^({pi, . . . + e, . . . \ ^^(P) and 

^Vhi^CP)' = ^^/«'^(^) \ ^/^^({Pl> . . . + e, . . . Now we want to prove 

ttiat tlie sum of tfie second term of ([20)) and the last term of ([^^ is null. 
First of all, recall that the sum in the second term of (l20l) can be limited 
to the cells in the neighborhood of the h-th cell, namely Vi G Ih with Ih — 
{j e {1, . . . , to} I G A/bz, {Ph, V)}- This fact implies that M ^ (AV;^ {V)+ U Al^f (P)-) 

AVi^{P)+ U A\4^(P)-. Moreover it can be easily seen that AV/^{'P)+ C 
AVj^lr)- and AV;f (V)- C AV;;^(P)+ Vi e h- Indeed, any segment s £ 
AV;^ (:P)+ is such that s e ^^({pi, . . . + e, . . . ,Pm}) and s ^ 1^^(7'), and, 
for any infinitesimal perturbation of p^, it is possible only if s e Vj^i'P) and 



^With A/gj-, (pht'P) we represent the set of neighbors of p/-^ in V. The neighboring property 
is given w.r.t. the proximity graph QdCP)^ that is the Delaunay graph associated to the 
Voronoi partition induced by V. 
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d 

lim 



dph^ 



Ph + ■ . ,Pm}), hence if s G AVf^(V) . Therefore we have 
/(||<Z-p,;||)<^(g) dq 



AV/^ (V) + 



fi\\q-p.\\)^{q)dq 



Av;^(P)- 



fi\\q~p,\\)<Piq) dq 



hm 



1 



"=^0 l|e|| yUiei^ ^K^CPV 

= hni 

e^o lie 



Ikll \JAV,^,(vr 



/ (Ik- Pill) 0(9) dq 
f{\\q-p^\\)q>{q)dq- 



AVi^^iV) + 



f{h-P^\\)<l>{q)dq 



f{\\q-p.\\)^{q)dq 



The conclusion follows from the fact that lime^o {^Vi^{P)+ U AV/^^iVy) = 
dV^{V) and Vg G dVf^ {V) \\q — Pi\\ — \\q ~ Ph\\- As concerns the second-last 
term of recalling again that lim._,o {AV^^{V)+ U AV^^{V)-) = dV^ {V) 
and the Assumptions 1161 we can write 



lim -—j- 
< lim — 



AV^^(V) 



lim 



dx 



[0,diam(AA)] 



{n\\q-{ph+e)\\)-f{\\q-ph\\))cjy{q) dq 
df 

M ll'/'ll[0,diam 

[0,diam(AA)] 

[0,diam(^)] {^ViiiP)^) ^ 0, 



dx 



dq 



where diam(A/') = maxp^ggjv" \\q ~ p\\- The same argument holds for the term 
with AV^('P)~, hence we have 



{fih-iph 



dph ' ||e||jyAr(p) 

-^!{\q-Ph\)^{q)dq 
V,^{V)OPh 

fc=i 



f{\\q-ph\\))cj,[q) dq 



with A/ft (P) the number of segments in V;^(P) and Ifc — j f {\\q — p^W) 4> [q] dq 

J St 



flfc 



and Sfc = [afc,6fc] G V^('P). If we choose the parameterization 7fc(i) 
(&fc — aA:)i, ^ G [0, 1] for Sfc we can apply the Theorem [T^ in appendix. Recall 
that in this case v{x,q) = \\q ~ x\\, hence the equation ||7fe(i) — Ph\\ — Ra ~ 
may have at most two zeros at t'^ ^ and 2 ^ {!'■•■ ' I* is worth noting 
that assumptions iii) and iv) play here the same role of assumptions i) and 
ii) in Theorem [TOl Therefore, from the definition of /(•), equation ([T]), we have 
the thesis. ■ 

In order to define a gradient-like algorithm, also in this case, we must relax 
Assumptions [T51 First of all, focus on the orthogonality assumption. It has been 
introduced to avoid the presence of entire segments in the boundary of a cell. 
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because these configurations induce problems in the definition of the gradient 
(they represent points on which the gradient may assume different values). Even 
in this case we opt to use the lexicographic rule in order to univocally assign a 
segment on the boundary to only one cell, and, again, we consider a discrete- 
time dynamics for the gradient-like algorithm. 

Using the lexicographic rule, we re-define the Voronoi cell as follows 

Vt^{V) = {g e I ||g < Ik - P.- II ^Pj A 

\\q — Pi\\ < \\q — PjW if j < i w.r.t. the L.O.} 

and verify that the expression ([TS]) for H{P) is still formally correct. We remove 
the orthogonality hypothesis by adding to (IT51) an Ik term for each segment 
entirely included in the boundary of a Voronoi cell. This fact does not change 
the expression (jl9p . since, with the new definition V^^(P), Mh{V) accounts now 
also for segments on the boundary. 

The relaxation of the other assumptions would imply some discontinuities in 
the integration domain induced by the discontinuities of the function /. These 
discontinuities, without additional assumptions, would prevent us from guaran- 
teeing H('P) to be monotonically nondecreasing along the evolution of V given 
by the gradient dynamics. Hence, we assume now / to be continuous and piece- 
wise differentiable. Being / continuous, the second term in Ik in (|19p is null. 

As made in the previous section, the directional derivative must be univo- 
cally defined on the vertices. To this aim, we use the expression given in 
Definition [TH but with reference to the formula for the directional deriva- 
tive in a point in the interior of a segment. Using these definitions we can state 
the following theorem. 

Theorem 18 Consider the following discrete-time evolution for the sensors' 
positions 

p{k+i) = p{k) ^ SkDHiV^''^), (23) 

where the h-th component of DTi is given by U6]) and Dwij'H{'P)[ph\ by il9\) 
and % : A/"™ IS^ as in // /(•) has locally bounded second derivatives, 

then, for suitable 5k, V'^^^ lies in a bounded set and 

i) 'H{'P'^^^) is monotonically nondecreasing; 

ii) ■p'-'"'^ converges to the set of critical points ofH. 

Proof. As long as sensors' configurations not violating orthogonality assump- 
tion are considered, the gradient is smooth and the proof is canonical. In the 
case of discontinuity points, segments belonging to the boundary of a cell can 
change allocation during sensors' motion. Hence, we can proceed as in the proof 
of Theorem 1101 and II 5 1 replacing barycenters with segments. In particular, being 
Sh a segment changing allocation, equations (jl3l) are now replaced by 

H,,(7'('=+1))=H«(P('=+1))-X^, 
with l[. = / -7^f{\\q — pi\\)4>{q)dq. Again, due to the fact that any point 



q ^ Sh is such that 



(fe+i) 
■Pi 



< 



(fc+i) 



and the monotonicity of /, we 



have I,* > , whereby the thesis follows as in the proof of Theorem (TU] 
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5 
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3.5 


9 


12.5 


13.5 


15 


14.5 


17 


20 
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2.3 


4 


5.5 


5 


8.5 


8.5 


7.2 


6.2 


10.5 


9 


7 




1.5 


1.5 


1.5 


2 


4 


1.5 


1.5 


1.5 


2 


2 


4 


ay 


1.5 


1.5 


1.5 


2 


4 


1.5 


1.5 


1.5 


2 


2 


2 



Table 1: Parameter values for the 11 Gaussian functions G{a,Cx,Cy,ax,<^y) = 



5 A case study 

In this section we apply the proposed two-step optimization procedure to a net- 
work representing a wing of the Amsterdam Schiphol airport. The network is 
made up of 63 vertices and 87 segments and the density function is the sum of 



with parameters assuming the values given in table [T] 

The network and a contour plot of the density function are shown in fig. [1] 
(darker colors denote preferential areas) . For sake of providing a clear graphical 
representation, the density function shown here is defined on M^, but the one 
used in all simulations is restricted to the network. 




Figure 1: Network representing a wing of the Amsterdam Schiphol airport and 
contour plot of the density function (j) used in simulations (darker colors denote 
preferential areas). 

The 50 sensors to be deployed have performance function f(x) = ^ ~ tanh 

where i? is a parameter considered as variable in the first step and as fixed in 
the second step of the optimization. Even if the previous function describes sen- 
sors with an infinite sensing radius, they will be represented as shaded circles 
of radius |i? to emphasize that the performance function assumes values lesser 
than 0.01 for larger distances. 




making up the density function (j). 



11 Gaussian functions of the form G( d , Cx ) Cy , (Ta; , (Ty ) — CI 6Xp 
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5.1 First step 



The first step of tlie optimization is performed on a collapsed network with 
collapsing factor r = 0.3 (see the small dots along the grey network in fig. ^ 
a,-c)). Sensors are grouped in 10 clusters of 5 elements each, and each cluster 
is represented as a single sensor. Clusters set initially i? = 10 and decrease 
linearly its value up to 1 during the simulation (compare fig.[5]-a) with fig.[21-c)). 
As apparent by the flows in fig. [2]-b), clusters are allowed to move in M?. 

It is important to recall that, both the variation of the sensing radius and the 
unconstrained motion of sensors are allowed in the first step as it is performed 
off-line. This step makes use of the algorithm described in section 13.11 and is 
thought to provide a good starting point for the second step. However, if sensors 
are initially located on the network as in fig.[21-a), they can execute the first step 
independently, using partial or rough information of the environment, without 
moving, and then plan a route on the network to reach the previously computed 
final positions. Since final positions can be not on the network (see fig. [ll-c)), 
they must be projected on it to be reachable. Anyway, this projection has to 
be performed before the second step to provide a valid starting point. 

5.2 Second step 

The second step considers a full network with sensors having fixed radius R = 
1 and initially deployed in the positions shown in fig. ^a). Such positions 
are obtained by spreading randomly 5 sensors close to each cluster center and 
projecting them on the closest segments of the network. Sensors now can take 
real measures from the environment and perform the optimization on-line, while 
moving, according to the algorithm described in sectional They are constrained 
to move on the network as shown in fig.[31-b). Final positions (see fig.[31-c)) show 
how sensors, originally clustered, diffused to better cover preferential areas (see 
fig©. 

6 Conclusions and Future Works 

This paper focused on the problem of optimally deploying sensors in an envi- 
ronment modeled as a network. An optimization problem for the allocation 
of omnidirectional sensors with potentially limited sensing radius has been for- 
mulated. A novel two-step optimization procedure based on a discrete-time 
gradient ascent algorithm has been presented. In order the algorithm to not get 
stuck early in one of the many local minima, in the first (off-line) step, sensors 
are allowed to move in the plane. Moreover, a reduced model of the environ- 
ment, called collapsed network, and sensors' clustering, are used to speed up the 
first optimization. The positions found in the first step are then projected on 
the network and used in the second (on-line) finer optimization, where sensors 
are constrained to move only on the network. 

The proposed procedure can be used to solve both static and dynamic de- 
ployment problems and the first step alone can provide solutions to location- 
allocation problems involving facilities located in the interior of the network. 

A main future research direction will consider the integration of classical 
Operative Research methods with the present gradient algorithm. In particular 
the first optimization could be addressed by adapting methods and heuristics 
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developed for the solution of the multisource Weber problem The aim 

is to build an overall global optimization technique to solve location-allocation 
problems of large dimensions with many facilities. Moreover, future research, 
more related to deployment problems, will consider other sensor's models such 
as those with limited sensing cone. 



Theorem 19 Let (p : M.'^ x ]R_|_ M. be a smooth function w.r.t. its second ar- 
gument and a non-increasing, piecewise differ entiahle map with a finite number 
of bounded discontinuities at Ri, . . . , G Ri < ... < Rn , w.r.t. its first 
argument. Let : x — >■ R_|_ be a continuously differentiable map w.r.t. 
both its arguments. Let s ~ [a, b] C be a segment and assume that 



i) v{x,a),v{x,b) <^{Ri,...,Rn}; 

ii) Vq e [a,b], ifVv{x,q) ■ (6 - a) = ^ v{x,q) i 

then 



,Rn}, 



jip {v {x,q),q)dq 



[o,i]\{t 

N ki 

i=ij=i 



d , , , , , , , , dvlx, 7(t)) 



\b-a\\dt 



diy{x,j{t)) 



dx 



X — X 

t=tij 



\b-a\\ 



where 7(t) = a + (b — a)t, t G [0, 1] is a parameterization for s and ti 
[0, 1], j G {1, . . . , ki\ are the zeros of v{x,^{t)) — Ri — 0. 

Proof. By using the Dirac's delta formalism we have 



T" ( ^{^{x^q) ,q)dq 
dx J, 



= / Q^V'(^(^'l)^l)dq 



— ip{u{x,q),q) 



dv{x,q) 



u(t{Ri,...,RN} 



dx 



N 

Y,{p{Rt.q)~V{Rl,q))6{v^R,) 

J s i^ 



dq 
di^{x,q) 



dx 



dq. 



With the chosen parameterization of s, the equation u{x^^{t)) — Ri — may 
have ki zeros G [0, 1], j G {1, . . . , fc^}. Thanks to i) and ii) the set of zeros 
tij does not change cardinality for x near a;, thus, tij depends smoothly on x. 
Recalling that for every arc 7 : [l,u] —^V we have 



f{x)dx 



f{^{t))m)\\dt, 
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and that for the special case of s ||7(f)|| = ||6 — a||, the derivative becomes 

^ [(fiif {x, q) , q) dq 

"'^ Js x=X 

\\b-a\\dt 



i^M., .(*))..(*)) 



\\b — a|| dt, 



•^[0,1] i=l 

from which, using the property / f{x)5{x — x)dx = f{x) for x G \l,u\, we 
have the thesis. ■ 
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a) 




c) 

Figure 2: First step of the optimization procedure for the deployment over a 
coUapsed network of 50 sensors grouped in 10 clusters: a) initial positions of 
clusters (sensing radius R — 10) and related Voronoi partition; b) gradient 
ascent flows (clusters allowed to move in R^); c) final positions, not necessarily 
belonging to the network, (sensing radius R = I) and related Voronoi partition. 





c) 

Figure 3: Second step of the optimization procedure for the deployment of 
50 sensors with R — 1 over a full network: a) initial positions obtained by 
spreading randomly 5 sensors close to each cluster position found in the first 
step (see fig. ^c)) and then projecting them on the network, related Voronoi 
partition; b) gradient ascent flows (sensors constrained to move on the network); 
c) final positions and related Voronoi partition. 
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